Optical spectra of a quantum dot in a microcavity in the nonlinear regime 
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The optical emission spectrum of a quantum dot in strong coupling with the single mode of a 
microcavity is obtained in the nonlinear regime. We study how exciton-exciton interactions alter 
the emission spectrum of the system, bringing the linear Rabi doublet into a multiplet structure 
that is strongly dependent on the cavity-exciton energy detuning. We emphasise how nonlinearity 
can be used to evidence the genuine quantum nature of the coupling by producing satellites peaks 
of the Rabi doublet that originate from the quantized energy levels of the interactions. 



(N 






cd 



T3 

C 

o 
o 



> 

On 
On 

en 

o 

oo 

o 



X 



I. INTRODUCTION 



Since the first experimental realization of atomic cavity 
quantum electrodynamics (cQED)fi a wealth of new re- 
sults have been obtained in the physics of strong coupling 
of a zero dimensional quantum system to a single pho- 
ton modei^^i^ The problem is important not only be- 
cause of the fundamental aspects brought forward by the 
interaction of material systems with photons^ but also 
because of the potential application of cQED to quan- 
tum information processing!^ In semiconductors, strong 
coupling is routinely achieved in planar microcavities — 
where in-plane excitons couple to photons with matching 
momenta — since the phenomenon was experimentally ev- 
idenced in 1992 iiS Strong coupling in 2D now provides 
the basis for investigating exotic phases of matter in these 
structures such as Bose-Einstein condensate a 11 ' 12 ! 13 ! 14 or 
supcrfluids i 15 ' 16 ! 17 ! 18 One system that has recently at- 
tracted particular attention is the zero-dimensional ana- 
log that consists of a Quantum Dot (QD) coupled to a 
single microcavity modei 19 ' 20 i 21 ' 22 i 23 The material excita- 
tions in the QD are excitons, that is, bound electron-hole 
pairs. Owing to their spatial confinement and energy 
level discretization, they can be brought in strong cou- 
pling with the single mode of a microcavity, such as that 
offered by a pillar (etched planar cavity)^ the defect of 
a photonic crystal 21 - or the whispering gallery mode of 
a microdisk ) 22 ! 23 among others. In Refs. |20ll2ll|22| | , such 
structures have demonstrated the Rabi doublet in their 
optical spectra, which is characteristic of the mode an- 
ticrossing that marks the overcome of dissipation by the 
coherent exciton-photon interaction. 



Steady progress is being made towards a better quan- 
tum strong-coupling as well as towards its external con- 
trol. Passive photonic crystals with breathtaking values 
of Q = 2.5 x 10 6 have now been reported^! limited only 
by structural imperfections. As the fabrication process 
will improve, a gain of another order of magnitude in 
the quality factor is allowed by the available theoretical 
designs. As for the active element, which inclusion low- 
ers the quality factor, prowess have been achieved in the 



engineering of the hetcrostructurc allowing to detcrmin- 
istically position a QD inside a photonic crystal to within 
25nm accuracy, and thus place the dot at a maxima of 
the light intensity, along with an etching technique of the 
holes of the photonic crystal to match spectrally the QD 
and cavity mode emission^ Micropillars etched out of 
Bragg mirrors allow a fine control of the radius and even 
of the shape of the cross-section, that can be optimised 
to reach high values of the quality factor, with figures 
of 1.65 x 10 5 reported for radius size of 4^m<2& One ap- 
pealing feature of these structures is their ease of access 
both for the excitation and the emission, the latter be- 
ing perpendicular to the sample surface. The mode vol- 
ume V in these structures can be further reduced while 
maintening the quality factor in one polarization by us- 
ing elliptical pillars, and record values of Q/VV can be 
achieved in this way— At the same time, the density of 
self- assembled QDs in the active medium has been suc- 
cessfully reduced over the years, with figures of 10 10 cm~ 2 , 
and the possibility to grow large dots (with lens shape 
of w 30nm) so as to provide a large oscillator strength. 
With microdisk resonators, the whispering gallery modes 
of a thin disk supported by a column provide the high Q 
modes of the cavity. 

These milestones open the way to new research in- 
vestigating the "true" quantum nature of the exciton- 
photon coupling. 2 ^ Although the system exhibits strong 
coupling, it is not known in which quantum state it is 
actually realized. To be useful for quantum information, 
the system must deal with few quanta of excitation and 
display strong alterations of its characteristic for small 
changes in particle number populations. If the system 
scales linearly with the number of particles, it is essen- 
tially classical^ and adding or removing a single par- 
ticle will not change its behaviour. Therefore, probing 
the nonlinear response is desirable. Recently, in an at- 
tempt to validate the single-particle character of the Rabi 
doublet in a QD-microcavity system, an unexpected ad- 
ditional peak was observed that was attributed to renor- 
malisation of the exciton energy by charging centres close 
to the QDi2£ The observed antibunching allowed the au- 
thors of Ref . [30| to confirm the quantum character of this 



coupling, but much is left to understand of the nonlinear 
effects and the nature of the additional peaks, observed 
or predicted. 

In this text, we compute the spectra of emission for 
a pumping high enough to enter the nonlinear regime. 
With QDs in microcavities, two types of strong nonlincar- 
ities are expected^ both associated to the active mate- 
rial, i.e., the excitons. The first one comes from Coulomb 
repulsion of the charged particles, and is the one inves- 
tigated here, in the case where it is comparable to the 
coupling strength. The other comes from Pauli exclu- 
sion, that arises from the fermionic character of the un- 
derlying particles^ In Refs. 3lH33| . the effect of Pauli 
exclusion alone was investigated (Coulomb interactions 
were added in an as yet unpublished work.—) There it 
was shown that Pauli exclusion resulted in rich multi- 
plet structures depending on the size of the dot, with 
a Mollow triplet-like structure for small dot with strong 
Pauli repulsion, and a Rabi doublet for large dots. In 
the present work, we show that even for bosonic exci- 
tons, a multiplet structure arises in the emission spec- 
trum when intra-dot exciton-exciton interaction is ac- 
counted for. For the large dots that the model aims 
to describe, Pauli exclusion can be taken into account 
phenomenologically by a phase-space filling effect that 
screens the exciton-photon interactioni 35 i 36 i 37 The latter 
results in loss of strong-coupling and we therefore focus 
on the intermediate regime where such renormalization 
of the coupling strength can be neglected. We also ne- 
glect the spin-degree of freedom, in particular the sign- 
dependent interaction between same and opposite-spins 
excitons, respectively^ This allows us to focus on non- 
linear deviations and to neglect more complicated cor- 
relations effects of the multi-excitons complexes such as 
formations of bound pairs or molecules that would give 
rise to bipolaritons i 39 i 40 ' 41 Experimentally, this could be 
realized by using a circularly polarised pump. Strong 
qualitative changes are still observed, with a basic struc- 
ture that consists of a blueshifted peak on top of a dis- 
torted Rabi doublet. 

The paper is organised as follow. In Sec. UH we in- 
troduce the Hamiltonian of the system describing the 
coupling between cavity photons and interacting QD ex- 
citons. From this Hamiltonian, we extract the energy 
level structure. This already contains the essential in- 
formation to understand qualitatively the spectra that 
we compute through more refined methods in the next 
section. In Sec. IIII1 we use the density matrix formal- 
ism and solve the master equation in the steady state to 
compute the spectra using the quantum regression theo- 
rem. This method, which suffers in principle no limita- 
tion or approximation (beyond assuming Markovian dy- 
namics), can only be used in practice for the low pumping 
regime, as the computational cost goes like ./V 8 with the 
method that we use here, where iV is the truncation order 
(the maximum number of photons accounted for in the 
model). We also propose an alternative formalism for the 
calculation of the emission spectra based on the Kcldysh 



Green function technique and compare the results ob- 
tained by the two methods. In Sec. US we discuss the 
overall results and conclude. The Appendix contains the 
mathematical details of the calculation of the emission 
spectra using the Keldysh Green function technique. 



II. MODEL HAMILTONIAN 

We describe light-matter interaction in a large QD in- 
cluding exciton-exciton interactions with the following 
Hamiltonian: 

H = J~] £]c) c 'J + ^r( c i c 2 + 4 c i) + -jr44 c 2C2- (1) 

.7=1,2 

The operator Cj (rf) is the annihilation (creation) op- 
erator for the cavity photon (j ~ 1) and QD exciton 
(j = 2). They all follow Bose algebra. The parame- 
ter Vr is the coupling strength and U accounts for the 
exciton-exciton interaction. In all the following, Vr is 
taken as the unit. An important parameter of the sys- 
tem is the detuning A = e\ — e 2 . The reference energy is 
that of the exciton, that is set to zero, £2 — 0. In cqn JT]), 
we neglected the polarization degree of freedom^ thus 
addressing the case of a symmetric QD under normal ex- 
citation, or excited by a circularly polarised light. 

To obtain an intuitive picture on the optical spectrum 
of such a system, we first study the allowed transitions 
between the possible energy levels. That is, we first ob- 
tain the eigenstates of the system with n excitations by 
diagonalizing the Hamiltonian {!]). Note that the total 
number of excitations is conserved. In what follows, we 
refer to the set of eigenstates of the total particle num- 
ber c\c\ + c 2 C2 with eigenvalue n as the nth manifold of 
the system. 

The structure of energy levels at resonance is sketched 
in Fig. (fT]) up to a maximum of two excitations (i.e., up to 
the second manifold) for three different cases. First, the 
bare levels corresponding to non-interacting and uncou- 
pled (or weakly-coupled) modes (Vr = 0, U = 0). Sec- 
ond, the eigenencrgies arising from the coupling (Vr 7^ 0, 
U = 0), and finally, the blueshifted lines that result from 
including the interactions (Vr ^ 0, U ^ 0). All the levels 
but those in the manifold n = 1, involving only one parti- 
cle, change with the interactions because of the excitonic 
part of their corresponding eigenvectors. In order to keep 
track of the excitonic character of each level, in Fig. |(2J), 
we plot the excitonic component of the eigenvectors of 
manifold n = 2 and their corresponding eigenencrgies as 
a function of the interaction U. We can see that, starting 
from a situation completely symmetric between the pho- 
tonic and excitonic fractions, the higher level gets more 
and more excitonic-like with U and blueshifts strongly. 
The other two energy levels are only slightly affected, as 
follows from their more photonic character. This charac- 
terisation of the levels, which also depends on the detun- 
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FIG. 1: (Colour online) Energy levels of the eigenstates of 
eqn {TJ up to the second manifold (two excitations), all at 
zero detuning (A = 0), left panel for weak (or no-) coupling 
(Vr. = 0), central and right panel in strong-coupling, with 
right panel also including interactions U varying on the x 
axis. The transitions between levels account for the spec- 
tral features. Red lines correspond to the vacuum-field Rabi 
doublet, that turns into a single Green line in weak coupling 
regime. Blue lines superimpose to the Rabi doublet when 
higher manifolds are probed. Without interactions, U — 0, 
these transitions do not appear in the spectra. Transitions 
from n = 2 to n = 1 in presence of interactions are plotted 
in Fig. |3] as a function of the detuning and U. New qualita- 
tive features appear thanks to the interactions. Black dashed 
lines are new transitions previously forbidden, although they 
remain weak. 



ing, plays an important role when identifying the spectral 
lines, as we show in the following sections. 

The energy levels of Fig. {TJ) are broadened by the fi- 
nite lifetime and by the incoherent pumping mechanism. 
This will be taken into account accurately in Section HTT1 
At this stage, an imaginary part is added to the energies 
of the bare states of the Hamiltonian (JTJ). This is justi- 
fied a posteriori by comparison with an exact numerical 
method. 

Once computed the new complex eigenvectors and 
eigenenergies, we obtain the amplitudes of probabilities 
of loosing an excitation from a given manifold to the 
neighbouring one counting one excitation less. We con- 
sider the amplitude for the processes of annihilation of a 
photon, that is, the transitions corresponding to the nor- 
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FIG. 2: Excitonic component at resonance of the eigenvectors 
corresponding to each energy level of the manifold n = 2 [see 
Fig. (JTJ)] as a function of the interaction strenght U. Varying 
the detuning also changes the character of the lines. In inset 
is plotted the exact energies of the eigenstates as a function 
of U, that are sketched in the right-upper part of Fig. JT} . 



mal mode emission. If the coupling is strong enough, the 
optical spectrum S p h(w) can be approximated as a sum of 
Lorentzians corresponding to each allowed transition^ 
The emission lines are finally broadened with the sum 
of the imaginary parts of the eigenstates involved in the 
transition and weighted by its probability. 

These results are plotted for transitions between man- 
ifolds with 2 and 1 excitation as a function of the non- 
linearity strength U on Fig. ([3]) a and as a function of 
the detuning on Fig. (J3J)b. In the latter case, the exci- 
ton bare energy is kept constant, equal to zero, while the 
cavity mode is brought in or out of resonance with the 
exciton. This relative detuning between the bare exciton 
and photon is tunable experimentally through a variety 
of techniques that include changing temperature or ap- 
plying a magnetic field (shifting the energy of the dot 
with negligible perturbation on the cavity) to growing 
a thin film (shifting the cavity mode energy without af- 
fecting the dot). As this is a parameter easy to tunc 
experimentally, that brings many and specific changes in 
the system, a density plot of photolumincscence emission 
with detuning of the type of Fig. [3] is the kind of results 
we are aiming for. On such a plot, one can trace the 
intensity, width and location (that include the number) 
of spectral lines for smooth changes in the detuning. 

Peaks appearing in Fig. ([3])a correspond to the transi- 
tions plotted in Fig. (JTJ). They arc labelled in blue (colour 
online): lower lines are the transitions "A" and upper 
lines the transitions "B" . Comparing with the linear Rabi 
doublet, which is superimposed in red, we observe the 
aforementioned blueshift of both groups of lines. It is 
more important for the exciton-like mode (especially line 
B-2 at resonance when U 3> 1) while the photon- like 
mode has a better resolved fine-structure splitting. At 
various detunings (Fig. ©b), complicated structures 
are found with crossing or anticrossing of the lines, as 
shown on the figure. Lines with the same bare-excitation 




one photon as the final state. They appear dimly in 
Fig. ©. 

Thus, the manifold method allows an understanding 
of the composition of the optical spectra, as illustrated 
on Fig. $5§ where the spectral lines have been labelled 
according to their corresponding transitions on Fig. |T]). 
From this overall picture, the excitonic fraction is clearly 
associated to the blueshift. Which transitions enter the 
dynamical picture, and to what extent, is investigated by 
the exact dynamical treatment that we pursue in terms 
of the density matrix formalism. 

In the next section, to describe more realistically the 
experimental situation, we supplement eqn ([1]) with an 
external excitation and a decay mechanism for the QD 
exciton and cavity mode. This is done by supplying a 
reservoir of external photons coupled to the cavity mode 
(to account for the finite lifetime of cavity photons) and 
another one to the exciton mode (to account for the ex- 
citon decay into photonic modes other than the cavity). 
Both pumping and decay compensate and the system 
reaches a steady state. 



III. CALCULATION OF THE EMISSION 
SPECTRA 

The coupling of the system to the external world re- 
sults in a leakage of photons and of QD excitations, that 
imply that even when the initial state of the system is well 
known (pure state) , the system evolves into a mixture of 
states where only probabilistic information is available. 
A suitable treatment of such a dynamics involves the den- 
sity matrix operator p. Its temporal evolution is given 
by a master equation of the type: 



FIG. 3: (Colour online) Cavity emission spectra from second 
to first manifold transitions: (a) Resonant case as a funtion of 
interactions U (y-axis). (b) Case of fixed interactions (17 = 2) 
as a funtion of detuning (A = e\ — €2, in x-axis). The non- 
interacting case U — 0, where only Rabi doublet arises, is also 
shown (red superimposed lines) as well as the bare cavity and 
excitonic lines (dashed green lines) . Lines are labelled in blue 
corresponding to the transitions of Fig. (ffj. Parameters in 
both plots are: 62 = 0, Tx — 0.1, T2 = 0.01, all in units of Vr. 



(photon or exciton) character cross, whereas lines of a 
different character exhibit anticrossing. At large detun- 
ings, the bare photon and exciton modes (in green) are 
recovered but with an additional blueshifted bare exciton 
line. 

Satellite peaks arise at very low and high energies from 
transitions that are forbidden in the linear regime. They 
enter the dynamics through nonlinear channels opened by 
the interactions. The dashed arrows in the right panel of 
Fig. |T]) represent these two transitions, with two excitons 
as the initial state that release one excitation and leave 



d P ■[ 771 , ^r 

- =l [p,H} + -C p 



(2) 



with H the Hamiltonian dynamics already introduced 
in eqn |T]), and with the incoherent contributions in the 
form of Lindblad terms Co, defined as 



Cop = 20pO^ - O^Op - pO^O 



(3) 



that corresponds to the general de-excitation operator O 
and its effective decaying rate x- Explicitly, the escape 
of the cavity photons is accounted for by a Lindblad 
term C ci p, with a rate x = Fi. This parameter is in- 
versely proportional to the cavity quality factor Q: 



T 1 = e 1 /Q. 



(4) 



The spontaneous decay of the QD excited states into 
any other mode than the one of the cavity, as well as 
the non-radiative decay, are taken into account with the 
term C C2 and its associated rate T2. Such a rate is typ- 
ically much smaller than the cavity emission rate, IV 
Finally, the cavity or the excitons can be pumped with 



a continuous incoherent pumping. The associated Lind- 
blad terms now involve the creation operators. The mas- 
ter equation in our system therefore reads: 

— =i[p,H\ +—C Cl p+—C C2 p+—C c ^p+—C c ^p. (5) 

We have included two possible kinds of incoherent pump- 
ing, namely, the cavity pumping at rate P\ and the elec- 
tronic pumping (via the exciton) at rate P 2 - Experimen- 
tally, the most common practise is electronic pumping, 
as the photoluminescence from the cavity is usually ob- 
served and pumping it would hinder the measurement. 
P\ could be taken into account to describe an effective 
cavity pumping due to, for instance, other dots that are 
weakly coupled to the cavity modes and populate the cav- 
ity with photons^. In the case of direct cavity pumping, 
it would be more relevant to observe the exciton emission 
for photoluminescence, or to turn to coherent excitation 
and probe the transmission, reflexion and/or absorption 
of the cavity. Here we compute the cavity photolumi- 
nescence for the two kinds of pumping each considered 
on their own, for the purpose of comparison, keeping in 
mind that the more experimentally relevant case is that 
of electronic pumping. 

We first obtain the steady state p^ of the system, by 
solving dp/dt = 0£1 From the knowledge of this density 
matrix we can compute the mean value in the steady 
state of a general operator fi as ($7) = Tr{p' ss ^J7}. More 
interestingly, the optical spectra of the system can also be 
obtained from p( ss "> and the master equation ([5]) thanks 
to the quantum regression theoremi 45 ' 46 We write the 
steady state equation in the matrix form: 







dp a 

dt 



2 , Ma/3,a'/3'Pa'/3' 

a'fl' 



(6) 



where the labels a and (3 index the whole Hilbert space, 
namely, in our case of two oscillators, a — {n, to} and 
= {p 7 q}. As a result, M a p !a /pi is a iV 4 x TV 4 ma- 
trix where N is the truncation of each oscillator's Hilbert 
space. In the computations, we have checked that the re- 
sults were independent of this truncation once it is taken 
large enough. The emission spectrum in the steady state, 
defined as 

1 f 00 

Sp h ( W ) = -Re / e^ T (c t 1 (0)c 1 (r))dT, (7) 



can then be obtained from Lax's hypothesis that the 
two-time correlators follow the same regression equations 
than the average: 

s ab (cj) = lM-T,( M + iul ^L>p>P ( Ap' B e<*) ( 8 ) 

a'f3' 

where A and B are the creation and destruction operators 



(«o 



of the transitions and p\. a i ri is defined as 



that is, in the case of cavity (normal) emission, 
where A = c[ and B = c\ : 



\ SS ) / _|_ -. (,ss) 

PA,n,m;p,q ~ V P ' ^Pn,m;p+l,q J 
-£*n,m;p,q \/V ®m,q *hi,p— 1 ■ 



(10) 

(11) 



The drawback of this method is the high computational 
cost involved in solving eqn © to obtain M and to sub- 
sequently invert it in eqn ([8]). Therefore, we present here 
the spectra for the low pump regime. This is not a se- 
rious limitation for the time being, however, as we are 
primarily interested in small deviations from the linear 
regime and to probe the first ladders of the quantized 
energy levels. We use pumping rates of 0.01, yielding av- 
erage number of excitations of the order (clci) « 0.1 with 
probability to have two excitons of the order of 0.01. For 
these figures, a truncation at the fourth manifold (N = 3) 
is enough to ensure convergence of the results. The other 
parameters are fixed to the following values, motivated 
by experiments: Vr = 1 provides the unit (experimen- 
tal figures are of the order of tens of /xeV), Ti = 0.1, 
T 2 = 0.01 and U = 2. 





0.7 




0.6 






fc 


0.5 


■£>, 




ti 


4 


3 




ri 




d 


0.3 


--J 






0.2 




0.1 




().() 






(9) 



-2024 
A=ei-e 2 

FIG. 4: (Colour online) Mean number of photons (dashed 
line) and excitons (solid line) as a function of the detuning 
between cavity and excitonic modes (A = ei — £2), for the 
case of cavity pumping (Pi = P, P2 = 0) in thin red (x5) 
and electronic pumping (Pi = 0, P2 = P) in thick black. 
Regardless of the detuning, an approximately constant and 
equal population is obtained for the cavity intensity, due to 
the balance between the effective coupling-strenght and the 
exciton-population. An asymmetry is observed with detuning 
due to the interactions that bring the exciton closer or further 
to resonance with the cavity mode. Parameters: £2 = 0, 
U = 2, Ti = 0.1, T 2 = 0.01, P = 0.01. 

Mean numbers of excitons and photons are plotted on 
Fig. (H} , for the two cases of cavity (only) and electronic 
(only) pumping. Close to resonance, A ss 0, both pump- 
ing yield approximately equal exciton and photon pop- 
ulations (note that the cavity pumping case has been 
magnified by a factor five). Detuning the modes re- 
sults in a collapse (cavity pumping) or increase (elec- 
tronic pumping) of the exciton population, as could be 
expected. Regardless of the kind of pumping, however, 
the cavity population is approximately constant. In the 




FIG. 5: Cavity emission spectrum as a function of detun- 
ing A. The U = case for photons, which corresponds to 
the Rabi splitting, is also shown in red solid lines as well as 
the uncoupled (Vk = 0, U = 0) cavity and excitonic lines 
in green dashed lines. As the mean number of excitons is 
very low (smaller than 1) the probability of having more than 
two excitons is very low, only the second to first and first to 
zero manifold transitions appear. For high positive detun- 
ing, Coulomb interactions generate an additional peak close 
to 62 + U(ri2 — 1). Note that peaks originated from transitions 
between manifolds with n > 2 appear with a very small in- 
tensity, due to the non-zero probability to have three excitons 
in the system. All energies are in units of Vr. Parameters: 
e 2 = 0, U = 2, Ti = 0.1, F 2 = 0.01, P = 0.01. 



cavity pumping case, this is because the exciton gets de- 
coupled and thus the cavity is pumped at a constant rate 
(one can actually see a small increase in its population). 
In the electronic pumping case, this is because although 
the coupling decreases, the exciton population increases 
in proportion so as to feed the cavity with a constant 



flux of photons. In both cases, an asymmetry is notable 
with detuning, because the interactions bring the cavity 
and the exciton modes closer or further from resonance, 
respectively, coupling them more efficiently for positive 
detuning and therefore allowing a larger production of 
excitons in that case. As a result, the nonlinear branches 
of the actual spectra (i.e., branches others than the Rabi 
doublet) for positive and negative detunings, shown on 
Fig. §5§, are not exactly as those shown in Fig. Qb. The 
blueshifted peak is more clearly seen in the positive de- 
tuning case thanks to this exciton population asymmetry 
with detuning. However, an excellent qualitative agree- 
ment is obtained with the manifold method, if one super- 
impose the vacuum Rabi doublet to the lines arising from 
higher manifolds. Depending on the pumping scheme — 
cavity (a) or exciton (b) — only quantitative features are 
changed that consist mainly in different lincwidths and 
intensities of the branches, that are otherwise well ac- 
counted for by the manifold method (see also Fig. \§§. 
While comparing Figs. [5ja) and (b), it should be borne 
in mind how the total population changes with detuning, 
as shown on Fig. [5] For this reason, panel (b) has a more 
complex structure, but this is due to the higher manifolds 
that can be reached with the electronic pumping. It is 
in fact possible to identify the third to second manifold 
contribution by extracting the lines in Fig. ([5]) that do 
not appear in Fig. Qb. These lines are clearly weaker 
due to the very low (but not vanishing) probability to 
have three excitons in the system. The transitions from 
even higher manifolds are too improbable to be seen in 
the spectra for the pumping considered here. The main 
differences between the two pumping schemes, if equal 
populations can be considered by adjusting the pumping, 
are therefore to be found in the linewidth and intensities 
of the lines. The positions of these lines embeds the most 
precious indications on the physical system. 

On Fig. ©, spectra are displayed for particular de- 
tunings, in solid black for the case of electronic pumping, 
and thin red for cavity pumping. The electronic pump- 
ing, which we repeat is the most relevant case experimen- 
tally, yields the most interesting spectral shape. On top 
of the Rabi doublet, the interactions produce additional 
peaks that, at large detunings, arc clearly associated to 
the exciton (panels (a) and (e)). Three peaks, El, E2 
and E3 are obtained that correspond to one, two and 
three excitons coupling to the cavity mode, respectively. 
As these are brought in resonance with the cavity mode, 
the linear Rabi doublet dominates (essentially because 
the efficient coherent coupling collapses the exciton pop- 
ulation) and satellite peaks are observed, that betray the 
quantum nature of the system, as the emission originates 
from transitions between quantised manifolds. In the ab- 
sence of interactions, the Rabi doublet is always observed 
independently of the total number of particles (this is the 
definition of a linear behaviour). Therefore, interactions 
are useful to evidence a quantum behaviour linked to 
quantized energy transfers, in the spirit of such experi- 
ments as those used with atoms in cavity to demonstrate 
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FIG. 6: (Colour online) Spectra for different detunings corre- 
sponding to vertical "cuts" in Fig (JS|. Both cases, of cavity 
pumping (Pi = P, Pi — 0) and electronic pumping (Pi = 0, 
P2 = P), are represented in thin red and thick black respec- 
tively. At very large detunings ((a) and (e)), multiple excitons 
occupancy is observed through the peaks El, E2, E3. Close 
to resonance, these result in satellites surrounding the linear 
Rabi doublet, that dominates because populations collapse at 
resonance. Parameters: 62 = 0, U — 2, Ti = 0.1, T2 — 0.01, 
P = 0.01. 



quantization of the light field. 47 Here, nonlinear features 
are observed directly in the optical spectrum, whereas in 
Ref. [47| , time-resolved measurements were used to probe 
anharmonic oscillations of the Rabi flops. This represents 
a notable experimental advantage, as measurements with 
cw incoherent pumping arc typically easier to perform 
than time-resolved spectroscopy. 

We conclude this Section by discussing the computa- 
tion of the emission spectra in terms of nonequilibrium 
Green functions, that could be used when the size of 
the Hilbcrt space forbids numerical computations with 
a master equation. Keldysh Green functions are now 
routinely used in description of electronic transport in 
mesoscopic systems 4 ^ and we believe that their applica- 
tion for the description of optical properties of QDs in 
the nonlinear regime is of methodological interest. Be- 
sides, this approach allows the discrimination between 
mean-field and correlation effects, which is important for 
the understanding of the physical origin of the additional 
peaks in the spectra of photoemission considered above. 
It can be shown (see Appendix for details) that within 
a Hartree-Fock approximation the emission spectra con- 
tains of a Rabi doublet and no satellite peaks appear. 
This corresponds rather well to the case of the planar 
microcavities with embedded quantum wells where the 
number of excitations is large, but does not give an ade- 
quate description of the spectra of QDs where the total 
number of excitations in the system is small and corre- 
lation effects play a substantial role. Taking them into 
account within the Keldysh Green function formalism al- 
lows to reproduce the multiplet structure of the spectra, 
as shown in Fig. (J7]) for the same parameters and oc- 
cupation numbers as in the case of cavity pumping in 
Fig. (|4|) and (J5]). The Green function results agree with 
those obtained with the density matrix formalism, the 
main difference being an anticrossing between the satel- 
lite branches clearly seen at Fig. ([7]), which is probably 
related to the truncation scheme we use. 



IV. CONCLUSIONS 

We have analysed the optical spectra of a QD in strong 
coupling with the single mode of a microcavity at the on- 
set of the nonlinear regime obtained through an incoher- 
ent pumping of either the cavity directly, or of the exci- 
ton, with nonvanishing probabilities of more than one ex- 
citation. We have considered a large quantum dot dom- 
inated by exciton-exciton interactions. Exact dynamical 
results were obtained with a density matrix master equa- 
tion and the quantum regression theorem. The Green 
function technique was also applied to support the results 
and suggest possible extension to arbitrary high number 
of excitations. 

We showed that a deep physical understanding can 
be attained by identifying the spectral lines with the 
transitions between quantised energy levels of the non- 
linear Hamiltonian. The spectral shapes observed, that 
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FIG. 7: Cavity spectral emission as a funtion of detuning 
A obtained via Green function technique. Parameters are 
similar to those in Fig. ([5]). 



probe the first manifolds of multiple excitations, are ac- 
counted mainly by the Coulomb energy shift on top 
of the linear Rabi coupling, with crossings and anti- 
crossings of the lines with detuning, depending on their 
opposed (photon and exciton) or identical (photon or 
exciton) character. The quantitative behaviour of the 
linewidth and intensities of these lines depends mainly 
on the dynamics of the system, such as the populations 
of the bare cavity photons and bare excitons. Electronic 
pumping of the exciton — the case most commonly real- 
ized experimentally as far as external incoherent pump- 
ing is concerned — allows to achieve high populations at 
nonzero dctunings thanks to the small exciton lifetime 
and the reduced effective coupling. The balance of these 
two factors nevertheless retains an approximately con- 
stant cavity population. As a consequence, the optimum 
experimental configuration to observe nonlinear effects 
in the photoluminescence spectra is at an intermediate, 
nonzero detuning, for instance when the detuning is ap- 
proximately equal to the coupling strength. There is an 
asymmetry with the sign of the detuning due to the in- 
teractions that further helps in characterising the nature 
of the nonlinearity (e.g., to which extent it comes from 
the exciton-exciton Coulomb interaction). The presence 
of satellites with detuning demonstrates emission from 
quantized manifolds, and as such is a signature of the 
quantum regime. The spectral drift of these lines with 
detuning is a useful tool to explicit the exact form of the 
Hamiltonian that accounts for the exciton nonlincarities. 
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APPENDIX A: OPTICAL SPECTRA WITH THE 
KELDYSH GREEN FUNCTION TECHNIQUE 

In this Appendix, we present the details of the calcula- 
tions of the photoemission spectrum with Keldysh Green 
function technique. 

To model the pumping and decay of cavity photons and 
QD excitons within the Keldysh Green function formal- 
ism, one needs to introduce the Hamiltonian of the cou- 
pling of the photonic and excitonic fields to the reservoirs 
of in-coming and out-going external photons and exci- 
tonic reservoir. The latter corresponds to the pump of 
the exciton state and its decay towards the leaky modes 
of the system. The complete Hamiltonian can be thus 
represented as a sum of that of eqn |T]) and the reservoir 
Hamiltonian H cxt given by 

H cxt = ^(^2;k4;k rf 2;k + t2;kC 2 4 : k + f 2;k C 2 rf 2;k) 
k 

+ Yl ( w i;k rft i ;k T, d i;k»? + ti;k,r,cid\. kv + i 1;k ,,,cl d i;k7,) 

(Al) 



k 

rj— in, out 



The operator di^ v {d\. k ) annihilates (creates) one 
incoming (?y = in) or outgoing (r) = out) free photon 
outside the cavity, c?2:k annihilates a photon in a leaky 
mode, the parameters ti ; k,fj an d i 2 ;k are cavity-reservoir 
and exciton-leaky mode coupling strengths respectively 
and wi^k are the dispersions of the external photons 
and leaky modes. In the Keldysh formalism, the photon 
spectral emission [counterpart of eqn |7|)] is recast as 



S ph (Lu) = -RciG< 1 (cu) 

7T 



(A2) 



where G^oj) is the lesser Green function in the fre- 
quency domain for the photons inside the cavity. In 
order to calculate the correlator in eqn (|A2|) . we de- 
fine the retarded Green function GJj (£,£') = — i9(t — 

ii)([cj(i);4(t')])^ The equation of motion of G r 3l {t,t') 
follows directly from the Hciscnbcrg equation for the op- 
erators Cj(t) = i[H + H cxt , Cj(i)}: 

dm ~ ^)G r 3l (t,t<) = 5(t - if)S il + V R G^(t,t<) 

+ 6 j2 UG*™( t ,t) + J dtiE5,(t > ti)Gt,(ti,0 , (A3) 



where 



GtfSftO = -iOit-mclWc^c^cjit')}) (A4) 



and 



^ = %Ei f i^i 2 3 1; ( S(^ti) 

kv; 



^EfekWftti), (A5) 



with 1 = 2, 2 = 1. The functions 

Sl£(*>*l) = -W(t-tl)<[dl;l B ,(t)d t 1;kl ,(t , )]> 

<$?(*> *0 = -W(*- *i)<[d2 ; k(*)4 S k(0]> (A6) 

correspond to the retarded Green functions for free par- 
ticles outside the cavity (ingoing and outgoing photons, 
and leaky modes). To close eqn (|A3|) . one needs the 

(2) 

expression for G- ■ - ( (i, t'), which involves the third or- 
der Green function. The full procedure leads to an in- 
finite set of equations, coupling Green functions of dif- 
ferent orders G < -™- ) (in a so-called Bogoliubov chain). To 
close this system, one needs to apply a truncation proce- 
dure, factorising Green functions of some given order into 
products of lower-order correlators. The simplest way to 
close the Bogoliubov chain is to apply the Hartree-Fock 
approximation, which consists in breaking the correlator 
G r ---\ (t, t') into products of two-operator Green function. 
More specifically, based on the Wicks theorem, one can 
apply the following approximations: 

([c\(t)c 2 (t)c 2 (t),ct(t')}) = 2(c\c 2 )([c 2 (t),ct(t')}) 
{[cl(t)c 2 {t) Cl {t),cl{t')]) = n 2 ([ Cl (i),ct (£')]> + 



(4ci><[c 2 (*),4(t')]) 

([n 2 (t)n 2 (t)c 2 (t),ci(t')}) = 2n 2 ([n 2 (t)c 2 (t),cl(t')}). 

However, in the problem considered, this mean-field ap- 
proximation appears to be too rough. Indeed, in this case 
the resulting spectra consist of a blue-shifted Rabi dou- 
blet without any additional satellites (result not shown) . 
To take correlation effects into account, in the follow- 
ing, we keep the correlators exact up to the second order. 
In this case, the equation of motion for G r - mn t (t, t') reads: 

[im + *i ~ £ » - £ ™] Gfi,l(*> O = *(* - *%».».« 

- v ^i/*'* / ) + VhGg^(t,0 + V R G r £li(t,t>) 
+6 n2 6 mn UG r j % (t, t') - 5 j2 UG r $ mnJ (f, t') 

+5rn2UG^l mnA (t,t') + S n2 UG r J ( il nn ^t,t'), (A7) 
where 

b jmn ,l(t) = ([c-(*)Cm(*)Cn(*)>4(*)]) (A8) 



and 



tr(3) 



({c](t)4 n (t)c n (t)c p (t)c q (t),c](t')}) 



According to the Wick's theorem for G r ^ 3 \ the follow- 
ing truncation procedure holds: 



G 



r(3) 



(t,t') = -id(t-t') X 



{(c](t)c n (t))([cl(t)c p (t)c q (t),4(t')}) + 

(c](t)c p (t))([cl(t)c n (t)c q (t),cl(t')}) + 
(c](t)c q (t))([cl(t)c n (t)c p (t),c](t')]) + 

(4n(t)c n (t)){{c](t) Cp (t)c q (t),4(t')]) + 

(cUt)c p (t)){[c](t)c n (t)c q (t),4(t')}) + 

(cl(t)c q (t))([c](t)Cn(t)c p (t),cl(t')])}. (A10) 



Substituting eqn (|A10[) in eqn (|A7[) we obtain a new 
equation for G r ^ only, 



,d_ 

dt 



-n c m 



G r ^ n ,(t,t') = ^-t%mn,l 



+<5 n2 <5„ m f/G;i 2 n ) ,(i,t') - <J,-aJ7 x 

{2niiGj^ fI (*,f ) + 2n im GTg ) jJ (t,f ) + 2n in G$>(«0} 

+<5 rn2 [/{2n jm G^ rU (i,t')+^G^ m ,(t,t') 
+ 2n mm G^ ]( (t, f ) + n mn G r /± A (t, t')} 

+8 n2 U{n jm G r ^l l {t,t') + 2n jn G r n % ntl (t,t) 

+ n nm G r ^li{t, + 2n nn G r W ntl (t, t')} , (All) 
with riji = (CjCi). Defining the Green function vector 

(-"•(2) _ (r r(2) r r(2) r r(2) ~r(2) ^r(2) ^r(2) ^ T 
±2-1 ~ V- 1 222. 1^122, 1 1^212, 1^121, 1 i^lll, 1 1^211,11 

(A12) 
(T is for transpose) , eqn (|A11|) can be written in a matrix 
form 

Gf\t,t') =Q r0 (t,t')b l + f dtiQ r0 (t,hM G r [ 2) (h,t'), 

(A13) 

where the vectors are underlined once and matrices twice. 
The matrix of zero order Green functions G^ -* is given 
by 

G/W(i,i') = mae(£ \£ \g? \g? \g? \gl™) 

(A14) 
where Diag(r) is the matrix with elements of r on its 
diagonal, zero elsewhere, and g n is defined according 
to [i§- t - e„] gn {0) (t, t') = 6(t - f). The matrix M is the 
sum of two terms, M = M + M , where 



M, 



(A9) 



/ -Vr 2V r \ 

-V R -A 2V R 

Vr A -Vr Vr 

Vr -Vr -A Vr 

2V R -Vr 

V 2V R -V R A ) 



(A15) 



10 



and 



and 



/ 1 + 6n 22 



Mu = u 



Qn 12 


nn 







1 



+ 6n 2 2 













2ni2 2n 2 2 



-4ra 2 i 



Finally, the vector b ; is defined as 



\ 





~2n 22 / 
(A16) 



25 2L n 22 

2$2ini2 
S 2 in 2 i + Snn 2 2 
S 2 inu + 5uni2 

2Sunii 

2Sun 2 i 



(A17) 



The quantities riji can be computed sclf-consistcntly 
together with the lesser Green function through the ex- 
pression riji = (i/2n) J GfAuj)duj. For simplicity, here 
we take riji as numerical parameters. That is, we sub- 
stitute the effect of the coupling to external reservoirs 
for an effective decay added as an imaginary part to 
the bare cavity energy ei — iTi/2 and to the bare ex- 
citon energy e 2 — iT 2 /2. The effect of the pump is par- 
tially included by adjusting nn and 7122 to the values 
calculated in Fig. (HJ. It is possible to show that the 
Keldysh-contour-ordered Green function G^„ i( T , T ') = 
—i([TccAT)c m (T)c n (T), c](t')]) has an equation formally 
similar to eqn (|A13|) 4£ Then, applying analytical contin- 
uation rules, we find in the frequency domain the equa- 
tions for the retarded (r) and the advanced (a) Green 
functions as well as for the lesser (<) Green function, 



£<<V) 



.<(0), 



+ g/^hm G a / 2) ( 



g r(0) (w)M G< (2) (w) 

a(2) 



(A19) 



In order to take into account the finite lifetimes of the 
exciton and photon modes, a phcnomenological broaden- 
ing was included in the lesser Green function of the free 
particles, which reads: 



[G<(°>H],=-2* 



2 



(u, 



e;) 2 + (£) 2 "' 



(A20) 



where j = 2 for the first three and j = 1 for the last three 
diagonal elements of the matrix G < ^°'(uj). 

Finally, we note that the equation for the retarded 
Green function of the first order can be written in a ma- 
trix form in the frequency domain, 



-V R 

1 
1 



u 



-v R 

€ 2 + i^ 







G 



r(2) 

222,1 



G\i 

Gr 
21 





r r(2) 
"222,2 



Gr 
12 
Gr 



(A21) 



while for the lesser Green function of the first order one 
has 



ijj - ei + i-^ 
-V R 



-v R 



62+*^ 



G< /~*< 
11 "-"IS 

"21 ^2: 



E< 




Ga ria 

11 u 12 

Ga ria 

21 "22 



U 



r <{2) r <{2) ) (A22) 
"222,1 "222,2 



G 



r,o(2) 



{uj) = G r > am {u)h 



G r ' a(0) (oj)M G r f (2 \uj), 



(A18) 



Solving the set of eqns (|A18[) - (|A22[) we obtain the emis- 
sion spectrum of the system, cf. eqn (JA2I) . 
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